arXiv:1503.00864vl [q-fin.PR] 3 Mar 2015 


Affine LIBOR models driven by real-valued affine 

processes 

Wolfgang Muller, Stefan Waldenberger 

Keywords: LIBOR rate models, forward price models, affine processes, 
volatility smile 

Abstract: The class of affine LIBOR models is appealing since it satisfies 
three central requirements of interest rate modeling. It is arbitrage-free, 
interest rates are nonnegative and caplet and swaption prices can be calcu¬ 
lated analytically. In order to guarantee nonnegative interest rates affine 
LIBOR models are driven by nonnegative affine processes, a restriction, 
which makes it hard to produce volatility smiles. We modify the affine 
LIBOR models in such a way that real-valued affine processes can be used 
without destroying the nonnegativity of interest rates. Numerical examples 
show that in this class of models pronounced volatility smiles are possible. 


1. Introduction 

Market models, the most famous example being the LIBOR market model, are very popular 
in the area of interest rate modeling. If these models generate nonnegative interest rates 
they usually do not give semi-analytic formulas for both basic interest rate derivatives, caps 
and swaptions. One exception is the class of affine LIBOR models proposed by Keller-Ressel 
et al. [14]. Using nonnegative affine processes as driving processes affine LIBOR models 
guarantee nonnegative forward interest rates and lead to semi-analytical formulas for caps 
and swaptions, so that calibration to interest rate market data is possible. 

This paper modihes the setup of Keller-Ressel et al. [14] to allow for not necessarily non¬ 
negative affine processes. This modification still leads to semi-analytical formulas for caps 
and swaptions and guarantees nonnegative forward interest rates, but allows for a wider class 
of driving affine processes and hence is more flexible in producing interest rate skews and 
smiles. Fonseca et al. [8] also propose a modification of affine LIBOR models. There driving 
processes are affine processes with values in the space of positive semidehnite matrices. The 
approach in this paper has the advantage that a flexible class of implied volatility surfaces 
can be produced with a much smaller number of parameters. 

The structure of this paper is as follows. In section 2 affine processes and their properties are 
reviewed. Section 3 introduces the necessary notation and market setup and reviews affine 
LIBOR models. It concludes with some comments on practical implementation. Section 4 is 
the main section of this paper. The first part presents the modified affine LIBOR model and 
semi-analytical pricing formulas for caps and swaptions are derived. The second part then 
gives some examples of usable affine processes with numerical calculations. 
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2. Affine processes 

Let X = {Xt)o<t<T be a homogeneous Markov process with values in D = ]R>q x M” realized 
on a measurable space (n,M) with filtration {lFt)o<t<T, with regards to which X is adapted. 
Denote by P^[-] and [•] the corresponding probability and expectation when Xq = x. X is 
said to be an affine process, if its characteristic function has the form 

E® [e“'^*] = exp ((/)i(u) + • aj), u£iR^,x£D, ( 1 ) 

where 4> ■. [0, T] x iM'^ —> C and V' : [0, T] x iM'^ —> with iM'^ = {« G : Re(M) = 0} and • 
denoting the scalar product in By homogeneity and the Markov property the conditional 
characteristic function satisfies 

E^ [e^'^^\Fs\ = exp {4>t-s{u) + 'tpt-s{u) ■ Xg). 

Accordingly affine processes can also be defined for inhomogeneous Markov processes (see 
Filipovic [7]), in which case the above equality reads 

E^ j-e^-M I ^ exp {(ps,tiu) + ' ^s), u £\R‘^,x G D, 

with : iM'^ —)• C and 'ips,t : —>■ C'^ for 0 < s < t. 

X is called an analytic affine process (see Keller-Ressel [11]), if X is stochastically continuous 
and the interior of the set^ 


V := 


f 


u E : sup E^ 

^Re(u)-Xs 

[ o<s<r 

I- 


< oo 


Vx E D 


( 2 ) 


contains 0^. In this case the functions cf) and ijj have continuous extensions to V, which are 
analytic in the interior, such that (1) holds for all u £ V. 

The class of affine processes includes Brownian motion and more generally all Levy processes. 
Since Levy processes have stationary independent increments, in this case = u and 

4>t{u) = tK{u), where k is the cumulant generating function of the Levy process. Ornstein- 
Uhlenbeck processes are further important examples of affine processes. They are discussed 
in section 4.2. 

The standard reference for affine processes is DufRe et al. [4]. There they give a characteri¬ 
zation of affine processes, where (p and ip are specified as solutions of a system of differential 
equations^. Of all the rich theory of affine processes the methods in this paper only use the 
specific form (1) of their moment generating function and the following property. 

Lemma 1; Let X be a one-dimensional analytic affine process and Re(M) < Re{w),u,w E V. 
Then Re{ipt{u)) < ipt{Re{w)), i.e. V'tjvnR is strictly increasing. 


can be described as the (convex) set, where the extended moment generating function of Xt is defined for 
all times t < T and all starting values x £ E. By Lemma 4.2 in Keller-Ressel and Mayerhofer [12] the set 
V is in fact equal to the seemingly smaller set |u E C"* ; 3® G int(Ii) : < oo| . 

^This also implies that X is conservative, i.e. ¥^{Xt £ D) = Lix £ D and 0 < t < T. 

®The fact that this characterization holds for all stochastically continuous affine processes was first shown in 
Keller-Ressel et al. [13] and later for affine processes with more general state spaces in Keller-Ressel et al. 
[!.')] and Cuchiero and Teichmann [3]. 
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Proof: The case D = M_|_ is already contained in Keller-Ressel et al. [14], In case D = M the 
lemma follows from the fact that by Proposition 3.3 in Keller-Ressel et al. [13] iptiu) = e^^u 
for some constant /3. □ 

Remark: If D = M_|_, it is known that both, ijj and are monotonically increasing (Keller- 
Ressel et al. [14]). With D = M this stays true for ■0, but not (/>, as the deterministic affine 
process Xt = xq — t shows. 

3. Interest rate market models 

Classical market models 

Consider a tenor structure 0 < Ti < • • • < Tjv < Tiv+i =: T and a market consisting of zero 
coupon bonds with maturities Ti,..., T^v+i. Their price processes {P{t,Tk))o<t<Tj, are as¬ 
sumed to be nonnegative semimartingales on a filtered probability space {lFt)o<t<T,l^), 

which satisfy P{Tk,Tk) = 1 almost surely. If there exists an equivalent probability measure 
such that the normalized bond price processes P{-,Tk)/P{-,T) are martingales^, the mar¬ 
ket is arbitrage-free. In this case we can define equivalent martingale measures for the 
numeraires P{t,Tk) instead of P{t,T) by 

_ 1 P{0,T) 

dQ^ Pin,T) PiO,Tk)’ 

In particular under the measure the forward bond price process P{-,Tk-i)/P{-,Tk) and 
the forward interest rate process 

are martingales. This is the basic market setup used throughout the rest of the paper. 

In the classical LIBOR market models forward interest rate processes are modeled as 
continuous exponential martingales under their respective martingale measure Hence 

forward interest rates are positive. Using driftless geometric Brownian motions as driving pro¬ 
cesses caplet prices are given by the Black formula (Black [1]) while swaption prices cannot 
be calculated analytically. Alternatively one can start with modeling the forward bond price 
processes P{-, Tk-i)/P{-,Tk) instead of forward interest rate processes. Using again exponen¬ 
tial martingales like a driftless Brownian motion it is then possible to analytically calculate 
caplet and swaption prices (see Eberlein and Ozkan [5]). The drawback of this approach is 
that forward interest rates will be negative with positive probability. 

Keller-Ressel et al. [14] proposed the affine LIBOR models, where forward interest rates 
are nonnegative while swaption and caplet prices can still be calculated semi-analytically, i.e. 
up to a numerical integration. The above approaches model the individual forward interest 
rate processes (resp. forward bond price process) with respect to the individual measure 


■^One can extend bond price processes to [0,r] by setting P{t,Tk) ~ p(t^t) ^ 

P{-,Tk)/P{-,T) is a martingale on [0,T] if and only if it is a martingale on [OjTfc]. Economically this 
can be interpreted as immediately investing the payoff of a zero coupon bond into the longest-running zero 
coupon bond. 
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under which they are a martingale. Contrary Keller-Ressel et al. [14] model the price 
processes P{-, Tk)/P{-,T), which are all martingales under the same probability measure 

Remark: Note that all models mentioned in this paper do not fully specify the whole term 
structure, but only part of it. In order to price derivatives not contained within the specified 
tenor structure it is necessary to specify some kind of interpolation scheme. Arbitrary in¬ 
terpolations may lead to arbitrage, however one can always choose an interpolation method, 
such that the model stays arbitrage-free (Werpachowski [17]). 


The affine LIBOR models 


This section presents a summary of the affine LIBOR model introduced in Keller-Ressel 
et al. [14]. On the filtered probability space {Ll,A,{Pt)o<t<T,Q^) consider a nonnegative 
analytic affine process X with a fixed starting value xq G 1^>o- Foi" tenor structure 
0 < Ti < • • • < Ttv < Tjv-i-i =: T define for fc = 1,..., and 0 <t <Tk 


P{ti 7fc) — \pUk-XT\ -p i. _ 4'T-t{uk)-\-i'T-t{uf:)-Xt 

P(t,T) •“ [e l-AiJ-e 


Uk > OjUfc G V, 


(5) 


where E*® [•] denotes the expectation with respect to a probability measure^ Q. These price 
processes are martingales and the resulting model is arbitrage-free. 

Writing 

-P(t, Tfc-i) _ P{t,Tk-i) j P{t,Tk) 

Pit,n) ~ Pit,T) / Pit,T) ^ > 


in (4) shows that forward interest rates being nonnegative is equivalent to normalized bond 
prices of (5) satisfying 


P{t,T,) ^ ^ P{t,Tp) ^ 

P{t,T) P{t,T) - • 


(7) 


T 

Since for x > 0, e“ ^ is monotonically increasing in every component of u the monotonicity 
for normalized bond prices in (7) is satisfied as long as ui > ... uw > 0. 

The parameters in (5) should be determined, so that the starting values of normalized 
bond prices P(0, Tfc)/P(0, T) = exp {(pTiuk) + V’T('Wfc) • xq) fit the initial term structure in¬ 
ferred from actual market data. For most affine processes every term structure can be fitted 
and for currently nonnegative forward interest rates this can be done using an decreasing 
sequence ui > • • • > up > 0 (see Keller-Ressel et al. [14]). 


Remark: Since X is nonnegative, the k-th normalized bond price is not only greater equal 
to one, but is bounded from below by the time-dependent constant exp {cpT-t{uk)), which is 
strictly greater than one. Accordingly in the affine LIBOR models forward interest rates are 
bounded from below by a strictly positive time-dependent constant. 


Affine LIBOR models lead to nonnegative forward interest rates. Additionally this specifica¬ 
tion is appealing because the density processes for changes of measures are again exponentially 


® Since Xq is fixed, contrary to to section 2 any dependence of probability measures on the starting valne of 
the Markov process X will be suppressed from now on. 
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affine in Xt, i.e. inserting (5) into (3) gives 

^4>T-TiAuk)+'‘pT-TiAuk)-XT,. 

dQ^ P{0,Tk) 

Moreover normalized bond prices and because of (6) also forward bond prices are of ex¬ 
ponential affine form. It follows that the moment generating function of the logartihm of 
normalized bond prices under QP is also of exponential affine form and that calculation of 
caplet prices is possible via a one-dimensional Fourier inversion. If the dimension of the driv¬ 
ing process is one, swaption prices can also be calculated via one-dimensional Fourier inversion 
(see Keller-Ressel et al. [14]). Hence this approach satishes both, nonnegative interest rates 
and analytical tractability of standard interest rate market instruments. If the dimension is 
larger than one, the exact price of swaptions can only be calculated via higher-dimensional 
integration, the dimension of which is the length of the underlying swap. Alternatively Grbac 
et al. [9] provide approximate formulas for swaptions. 

Practical application of the afhne LIBOR model 

Although this framework is elegant from a theoretical point of view, a practical implementa¬ 
tion faces several difficulties which shall be discussed here. 

First, calibration of interest rates and implied volatilities cannot be separated. The initial 
term structure can be fitted using the Uk, but the parameters Uk also have a strong impact 
on implied volatilities. This can be seen by looking at the forward bond price 

1 

^ = exp {(t)T-t{uk-l) - (t)T-t{uk) + (ipT-tiuk-i) - ipT-t{uk)) ■ , (8) 

which is the random variable responsible for the payoff of a caplet. The driving process X 
influences the distribution of this random variable through two different channels. First via 
the parameters of the driving process itself and second via the parameters Uk (depending on 
X and the initial interest rate term structure). Hence for changes in the yield curve different 
parameters are required to reproduce the same implied volatility surface. If A is a Levy 
process, then as mentioned in section 2 iptiu) = u and it follows that the distribution of (8) 
depends on the difference Uk+i — Uk, which in turn is related to the steepness of the initial yield 
curve®. Hence caplet implied volatilities are especially sensitive with regards the steepness of 
the initial yield curve. 

Second, interest rates and volatilities of this model depend on the final horizon T. Changing 
the horizon T while using the same affine process X will lead to different results and there 
is no general way of rescaling the parameters of X to negate such an effect. This is rather 
counterintuitive, since extending the horizon of a model should not change the results for 
quantities already included with the shorter horizon. 

Third, the types of possible volatility surfaces is rather constrained in the fully analytically 
tractable one-dimensional case. For example, we were only able to generate volatility skews^. 


®This is similar for most affine processes, but is best visible for Levy processes. 

^The smile example of Keller-Ressel et al. [14], figure 9.2, using an Ornstein-Uhlenbeck process seems to be 
numerically incorrect for strikes smaller than 0.4. With the mentioned initial yield curve the underlying 
interest rate is always larger than the strike, which corresponds to a zero implied volatility, destroying the 
displayed smile. 
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This might be resolved by using higher-dimensional nonnegative processes. However, in multi¬ 
dimensional affine LIBOR models swaptions can no longer be calculated efficiently by Fourier 
methods. On the other hand allowing arbitrary affine processes destroys the nonnegativity 
of forward interest rates, a central property of affine LIBOR models. We propose a modifi¬ 
cation, that preserves the nonnegativity of forward interest rates without the restriction to 
nonnegative affine processes. 


4. The modified affine LIBOR model 

On the filtered probability space (H,^, {lFt)o<t<T,Q^) consider an analytic one-dimensional 
affine process X with a fixed starting value xq, i.e. the set V defined in (2) contains 0 in the 
interior. For u € V with —u G V consider the martingales M“, 

M“ := [cosh(uXT)| Xt] = ^ g</.r-t(-0+bT-t(-0^t^ , (g) 

By the symmetry of the cosinus hyperbolicus = M““, hence one may restrict u to be 
nonnegative. For the given tenor structure 0 < Ti < • • • < Tn < F/v-i-i = T and the market 
setup of section 3 define the normalized bond prices for A; = 1,..., iV and t < Tfc as 

:= M“^ Uke{vGV:v> 0, -u G V}. 

With being a Q^-martingale the model is arbitrage-free. For every x G M the function 
u I—cosh(rix) is increasing in u G M>o and satisfies cosh(rtx) > 1 so that if 


Ul > U2 > ■ ■ ■ > UN > 0, 


equation (7) holds and forward interest rates 


F^{t) 


Afc V ) ’ 


0 < t < Tk_i. 


are nonnegative for all t. To fit initial market data one has to choose the sequence {uk) so 
that Mg'" = P{0,Tk)/P{0,T). The following lemma gives the condition for the affine process 
X under which a given initial term structure can be reproduced and shows that the Uk are 
uniquely determined. 

Lemma 2: If 

P(0, Ti)/P{0, T) < sup [cosh(uXT)| Pq] , 

tiSV:—nSV 

then the model ean fit any term structure of nonnegative forward interest rates. Additionally 
there exists a unique decreasing sequence ui > ■ ■ ■ > un , such that 


P(0,Tfc)/P(0,r) = E'®^ [cosh(ufcXT)| Po] = 


If forward interest rates are strictly positive, the sequence is strictly decreasing. 

Proof: m{u) = E*®^ [cosh(riAr)| Pq] is a continuous function which is strictly increasing for 
u > 0. By the assumption of the theorem there exists u > 0 with m{u) > P{0,Ti)/P{0,T). 
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Furthermore m(0) = 1, which proves the lemma. □ 

Remark: Generalizing this approach to a d-dimensional driving process is possible by setting 


M“ = E 


n 

U=i 


cosh 


u 




u = > 0 . 


In this case it is guaranteed that M“ > for u > w, which guarantees the nonnegativity 
of forward interest rates. However, the option pricing formulas in the following sections do 
not generalize. 


As in the affine LIBOR model for a monotonically decreasing sequence (uk) forward interest 
rates are not only nonnegative, but bounded below by strictly positive time-dependent con¬ 
stants (the bounds can be calculated numerically). This is not a big issue if these bounds are 
close to zero, but has to be checked during the calibration process. 

In the modihed affine LIBOR model the change of measure to the T^-forward measure 
is given by 


P{0,T) Afg 

dQ^ P{0,Tk) 


( 10 ) 


Here is a sum of exponentials of Xt, while in the affine LIBOR model the corresponding 
term is a single exponential. This means that contrary to the affine LIBOR model the process 
X is not an inhomogeneous affine process under and it is not possible to calculate the 
moment generating function of the logarithm of foward bond prices under . Nevertheless 
it is possible to get analytical formulas for the prices of caplets and swaptions. 


4.1. Option pricing 

The derivation of the pricing formulas for caplets and swaptions is based on a method hrst 
applied in Jamshidian [10]. First caplets are dealt with, swaptions follow afterwards^. Note 
that if Uk = Uk-i the corresponding forward interest rate always stays zero. To exclude 
such pathological examples assume that the sequence (uk) is strictly decreasing. In this 
section random variables are often viewed as functions of the value of the driving process X. 
Specifically consider the functions : M —>■ M, 

X h-). M“(x) := - ^g</'T-t(“)+bT-tWa; _j_ g0T-t(-“)+bT-t(-“)^^ _ (H) 

The time t value of martingale M“ in (9) is then = M^{Xt). In the rest of the paper 
will denote both, the function and the value of the stochastic processes, where the correct 
interpretation should be clear from context. 

The payoff of a caplet for the {k + 1)**^ forward rate F^~^^{Tk) with strike K is 


,k+i [F’^+\Tk)-K) 


1 


P{Tk,Tk+i, 


-K\ = 




- K 


® Actually caplet prices coincide with prices of swaptions with only one 
between those two derivatives is the payoff time. 


underlying period. 


The difference 


7 











AfEne LIBOR models with real-valued processes 


where K = 1 + A^+iK. Since this payoff has to be paid at time the price of the caplet 
and the corresponding floorlet is 


Cpl{t,Tk,Tk+uK) = P(t,rfc+i)E^''''+' 


Flt{t,n,n+i,K) = P(t,Tfc+i)E^'''“+' 



IFt 

Tt 


Since price processes are martingales, the put/call parity holds and prices of caplets follow 
from floorlets and vice versa. Because Fourier analysis is easier for floorlets, where the payoff 
is bounded, formulas are derived for floorlets. 

Since the moment generating function of ln(M^^/My^^^) is unknown, Fourier methods are 
not directly applicable. However, the function x i—)• (x) has a unique minimum 

and is monotonically increasing moving away from this minimum. Using this one can get 
rid of the positive part and use Fourier inversion to calculate the above expectations. The 
above mentioned monotonicity is very fortunate and follows from a close interplay between 
the monotonicity of the sequence (u^) and the function V’ with properties of the cosinus 
hyperbolicus. Details are laid out in the proof of the following lemma, which can be found in 
the appendix. 

Lemma 3: For i = 1,... ,n let uq > Ui > 0, where for at least one i uq > Ui. Let Ci > 0 be 
positive constants. Define a function <7 : M —?• M by 

. Mfifix) 



Then g has a unique maximum at some point ^ G M and and is strictly monotonically de¬ 
creasing to 0 on the left and right side of f,. 


For floorlet valuation this lemma is not directly applicable as Uk > Uk+i, which is the wrong 
inequality. However, there is only one summand and the lemma can be applied to the inverse 
My*+^(x)/M^*(x). It follows that has a unique minimum at some point 

f and is increasing to inhnity to the left and right. Hence it is possible to write 


K - 


Mlpfix) ' 


= K- 


Mlpfix) 

M^x 


I{ki < X < K2] 


(13) 


where ki and K 2 are two uniquely determined constants satisfying ki < ^ < K 2 - If = C = ^2 
the payoff is zero, which corresponds to > K. This happens if the forward interest 

rate is bounded from below by K, which only happens for very low strikes K. Inserting (13) 
into the price of a floorlet it follows by a change of measure that 


Flt(t, Tfc, Tfc+i, K) = P{t, Tk+i)E^' 


k+1 




Tt 


P{t, r)E^^ I {ki < Xt, < K2} I Pt] ■ (14) 
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is the sum of exponentials of the random variable Xt,.- The expectation in 
(14) is calculated under the measure where the conditional moment generating function 



exp {(kt-siz) + 'ipt-siz)Xs) 


is known for z G V. Hence the expectation in (14) can be calculated via Fourier inversion. 
The Fourier inversion formula for terms of the above form is stated in Lemma 4, the proof of 
which is given in the appendix. 


Lemma 4: Assume that the function / : M —>• M has the representation 


f{x) = ^ C'fce’''"* I{ki < X < K 2 } , lim /(x) = lim f{x) = 0 , 

x],K,i X'\K2 

k 


where the summation is over a finite index set and the Ck and Vk are real constants. Then 
for i? G V n M the Fourier inversion formula 

E[f{Xt)\ Xs] = 11^ Re [Mx,\xS^u + R)f{u - xR)) du 

holds, where f is the analytic Fouier transform given by 

/(^) = - E z 7 ^ 0, z -ivk. (15) 

XZ Vk — xz\ J 

To calculate the price of a floorlet in (14) apply Lemma 4 to fk+i{XTfi) with 

fk+i{x) ■■= (kMl^l+\x) - M“^"(x)) I{ki < X < K 2 } . (16) 


Its Fourier transform is 

fk+i{z) = ^ {{1 + Xk+iK)hl\,^^{-xz,Uk+i) - hll,^^{-xz,Uk)) 

with 


.-e 2(z + ibT-fiu))V ) 


+ e' 


2{z + fir-tiu)) 

(tr-ti-u) i’T-tj—u) ^^(^z+TpT-ti-u))i^2 _ p(z+bT-t(-«))Kl 

2{z + xpT-ti-u)) 


(17) 


(18) 


The case of swaptions is similar. Consider a swap which is part of the tenor structure. That 
is, consider 1 < a < fi < N and the according interest rate swap with forward swap rate 


SaA^) 


P{t,T^) - P{t,T^) 

ELa+l^kPit,Tky 


Xk — Tk — Tk—i- 
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The payoff of a put swaption on the above swap with strike K is then 


/3 

P{Ta, Tk)Ak (K - 

k = 0!. + l 


/ /3 N 

PiT^,T^) + K Y, AkP{T^,Tk)-l 

\ k=a-^l / 

" P Y - 1 ' 


Mr 


La k=a-\-l 


'Mil 


Since the function (x)/ (x) + Ylk=a+i KA^Mlp^ (x) /(x) is of the form of Lemma 

3, it has a unique maximum ^ and one can find constants® ki < ^ < K 2 such that after a 
change of measure the value of a put swaption is 


where 


PutSwaption(t, Tq, T^, LC) = P(t, r)E‘® 




Tt 


fl^ix) = (x) - M““(x) + KA.MAix) \ I{ki < x < K 2 } . (19) 

V k=a+l J 

Again this of the form in Lemma 4 and in this case 

k=OL-\-l 

where is dehned in (18). The pricing formulas are summarized in the following 

theorem. 

Theorem 5: Let P G VnM. In the modified affine LIBOR model priees of a forward interest 
rate put and a put swaption are 

Flt{t,Tk,Tk+uK) = Re (^Mxrjx,{R + iu)fYA-^R)) d«, (21) 

PutSwaption(t,ra,r^, AT) = J Re (^A4xy^|Xt(P + in)/^^(u - iP)^ du, ( 22 ) 

The Fourier transforms repectively f^ are given in (17) respectively (20) for R ^ 
{0, Uk, Uk+i} respectively R ^ {0, u^, ..., up}. 


®As in the floorlet case if ki = K 2 then the forward swap rate is always larger than the strike. Note that 
Sa.fiifi) can also be written as Sc,pit) = with Wk > 0 (see e.g. Brigo and Mercurio 

[2]). It follows that if forward interest rates are bounded below by positive constants the same will be 
true for forward swap rates. This bound is then at most an average of the corresponding forward interest 
rates bounds and is therefore of the same order of magnitude, which for a meaningful model will be small 
enough. 
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In order to calculate respectively one has to find the roots ki , K 2 of the functions 


9k (x) : 

9aA^) ■ 


= K - 


J- cr 




(x) 


J a 


(x) 


+ 


0 

E 

k=OL+l 


KA, 


M^^Jx) 

'M-Hx) 


- 1 . 


(23) 

(24) 


By Lemma 3 this amounts to finding the roots of a function which has a single optimum 
and is monotonic when moving away from this optimum. Numerical determination of the 
roots of such well-behaved one-dimensional functions poses no problem. Having determined 
those bounds valuation reduces to a one-dimensional integration of a function that is falling 
at least like 1/x^ (depending on the moment generating function of the affine process), so 
also numerical integration is feasible. Note that besides caps, floors and swaptions, options 
like digital options or Asset-or-Nothing options can be calculated in a similar manner. 


4.2. Examples 

The first part of this section looks at the benchmark case of a Brownian motion, where 
everything can also be calculated in closed form. Afterwards Ornstein-Uhlenbeck processes 
are discussed. The section concludes with examples of possible volatility surfaces. 


Brownian motion 

Choose Xf = Bf, a standard Brownian motion starting in 0. The conditional moment gener¬ 
ating function is 

Mb^IbM = E = exp (^uBt + y (T - t)^ 


2 

Hence this is an affine process with cj)t{u) = \t and ipt{u) = u. Consider the time 0 price 
of a floorlet as given in (14) with t = 0. Since M“(—x) = M“(x) one finds that in this case 
k ,2 = K and Ki = —K, where k is the unique positive root of (23) if < 0 and k = 0 

otherwise. By (14) the floorlet price Flt(0, T^, T^+i, AT) is 


p(o, r)E‘^'^ 


^Ke cosh(rifc+iHrfc) - e cosh(«feHTj,)^ 


< k} 


By the symmetry of a Brownian motion starting in 0 

E[cosh(zHt)I{|Ht| < k}] = E[e^^n{\Bt\ < ac}] = E I {|Ht| < k}] 
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where denotes the cumulative distribution function of a standard normal distributed random 
variable. Hence 

FltiO,Tk,Tk+uK) = KP{0,T)e^U^^ - Uk+iV^) - ^ - n^+i 

-P(0,T)e«JJ _ $( _ ^ . 

Slightly more complicated formulas exist when B is replaced with a Brownian motion with 
constant drift and volatility and a starting value different from 0. 

Swaptions can be treated the same way. Let k be the unique positive root of (24) if > 

0 and K = 0 otherwise. Then 


PutSwaption(0, Ta^Tjs, K) = P(0, T)e^P 2 



up 





- P(0,r)e““? 





up 



0 

+ P{0,T) Y, 

A:=a +1 





- Uk 



Ornstein-Uhlenbeck (OU) processes 

The OU process X generated by a Levy process L is defined as the unique strong solution of 
(see Sato [16], section 17) 


dXt = -XXt dt + dLi, Xq = xq. (25) 

Then Y) := e^^Xt = xq + e^^Lg ds. Using the key formula of Eberlein and Raible [6] it 
follows that 


E[e 


uXtl 


= E 


exp 


e~^^uYt) = exp 


e ^^uxn -I- 


As„.\ 1 „ 


where k{u) = ln(E[Li]) is the cumulant generating function of the Levy process L. Hence 
this process is affine with 

= e~^^u and (ptiu) = [ K{e~^^u)ds = \ [ du. (26) 

^0 ^ J ^ 

By Corollary 2.10 in Duffie et al. [4] every affine process with state space M is in fact an OU 
process. Hence in the context of affine processes defined on the real line OU processes are the 
right class to consider. For application it should be possible to calculate the integral in (26) 
analytically. Two examples where this is possible are presented below. 


Remark: If L is a martingale, the process in (25) is mean-reverting to zero, however shifting 
the mean to 6 is easily done by using Zt = 0 + Xt. Then dZt = X{9 — Zt) dt + dLt and 

E [e“^*] = exp + 6u{l - e“^*)) -b ijjtiu)Zo'j . 
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Hence Z is again affine with and (/>f(u) = (f)t{u) + Ou{l — Note that this 

process is then generated by the Levy process Lt = Lt + 6Xt, i.e. the original Levy process 
plus an additional drift of OX. 

The first example is the classical OU process generated by a Brownian motion aB, where 
K,{u) = This process is described by 

dXt =—XXtdt + crdBt, Xq = x. 


The integral in (26) is 


Mu) 



k{vu) 

V 


dt> 


4A 


(1 -e 


-2Xt\ 


(27) 


With Brownian motion describing the continuous part of Levy processes., for the second 
example we consider a pure jump process, namely a Double T-OU process. T-OU processes 
are generated by a compound Poisson process with jump intensity A/3 (A being the same as in 
(25)) and exponentially distributed jumps with expectation value a. The limit distribution 
of this process is a T-distribution, which gives the process its name. As the generating 
compound Poisson process is strictly increasing, the generated P-OU process is a subordinator 
and stays above 0. In order to find an OU process with values in M consider the difference of 
two independent compound P-OU processes with parameters a"*",/3"'“, a“,/3“ and set 

A"*" = A/3"*“,A“ = A/3“. Then L = L'^ — L~ is a compound Poisson process, where positive 
jumps with expected jump size ^ are arriving at rate A"*", while negative jumps with expected 
jump size ^ are arriving at rate A“. 

The cumulant generating function of a compound Poisson process with exponential jumps 
is which is defined for u < a. Hence for u G (—a“,a^) the moment generating function 
of the combined process L is 


E[e“-^i] = E 


e“^i 


E 




= exp 


)tt^ + (/3+a — fd a'^)u 
(a+ — u){a~ + u) 


Inserting this into (26) straightforward calculations show that the function cj) for the resulting 
OU process is given by 


Mu) 


M + M 
2 


+ 



In 

In 



— e ^^u){a +e 
(a+ — u){a~ + u) 

— e~^^u){a~ + m )\ 

— u)(a~ + e~^*u) J 



(28) 


It is also possible to combine the two approaches by considering an OU process generated 
by a Levy process which is the difference of two compound Poisson processes plus a Brownian 
motion, all of which are independent. The resulting (p then follows by adding up the two 
functions (27) and (28) and for this process V = {u G C : —a~ < Re{u) < a"*"}. By the 
previous remark it is also possible to shift this process by 9. Such OU processes are used in 
the following numerical examples. 
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Figure 1: Implied volatility skew of caplets generated by an OU process with 
parameters A = 0.02, = 12, = 10, Z?"'' = 50,/3“ = 5, a = 0.3,0 = 

0.5, X = 0.7 and T = 10. 

Volatility surfaces 

With the OU process of the previous section it is possible to generate volatility smiles as well 
as volatility skews. For illustration we consider a term structure with constant interest rates 
of 3.5%. The tenor structure and therefore the forward interest rates are based on half year 
intervals. Implied volatilities are then calculated for caplets with maturities over a 5-year 
period and strikes ranging from 0.02 to 0.07. Figure 1 shows a skewed volatility surface while 
figure 2 shows a very pronounced smile, both of which are generated by an OU process of the 
just introduced type. As mentioned in the previous chapters forward interest rates in this 
type of model will be bounded from below. The bounds in these examples are at 1% for the 
forward interest rate expiring after half a year and decrease to basically 0% for the forward 
interest rate which expires in 5 years. Hence they are well within reasonable boundaries. For 
completeness an example of at-the-money implied volatilities for swaptions with maturities 
and underlying swap rates ranging from 2 to 7 years is displayed in hgure 3. 


Conclusion 

Classical interest rate market models are not capable of simultaneously allowing for semi- 
analytical pricing formulas for caplets and swaptions and guaranteeing nonnegative forward 
interest rates. One exception are the affine LIBOR models presented in Keller-Ressel et al. 
[14]. This paper modifies their approach to also allow for driving processes which are not 
necessarily nonnegative. Caplet and swaption valuation is possible via one-dimensional nu¬ 
merical integration. This allows for a fast calculation of implied volatilities for these types 
of interest rate derivatives. With the additional flexibility of real-valned affine processes this 
type of model is capable of producing skewed implied volatility surfaces as well as implied 
volatility surfaces with pronounced smiles. 
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Figure 2: Implied volatility smile of caplets generated by an OU process 
with parameters A = 0.02, = 50, a~ = 5, /S'*" = 50, /3“ = 10, d = 0, 0 = 

0, X = 1 and T = 10. 



Figure 3: Swaption implied volatilites generated by an OU process with 
parameters A = 0.02, a"'' = 12, = 10,/S"'' = 50,/?“ = 5,(T = 0.3,0 = 

0.5, X = 0.7 and T = 10. 
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A. Proofs 

Proof of Lemma 3; For a function f{x) denote its even and odd part by 

fix) = ^ifix) + fi-x)), r{x) = ^(/(x) - fi-x)). 

Note that if / is monotonically increasing, the same is true for f°. Then (11) can be written 
as 

M^{x) = - ^Q't>T-t(u)+'>pT-t{u)x _|_ 

= e^T-tO)+rT-tO> + ^T-tiu)x) 


and 


Mfix) 

W¥) 


— g(</'T_t(«i)-<AT_t(“o)) + (bT_tK)-V’T_t(«o))a: 


cosh(0y_j(ui) + 'il;^_^{ui)x) 

cosh((()^_^(uo) + V’T_t('Uo)x)' 


(29) 


If Ui = uq, then (29) is constant and has no influence regarding monotonicity or maxima. 
Hence from now on assume uq > Ui for all i. The function g of equation (12) can be written 
as 


aix) = 

i=l 


.A,.^aiX cosh(Hi + bjx) 
cosh(ilo + h^x) ’ 


where for i = 0,... ,n 


— ifr-ti'^i) ~ 4’T-ti'^o)), Bi — cj)^_^{ui), 

Ui = [ipT-tiUi) - 'll^T-tfo)), h = ll^T-tiUi)- 

Since is monotonically increasing (see Lemma 1), also ■0° is monotonically increasing. With 
'ijj°{0) = 0 it follows that bi > 0 for all i. Furthermore note that a* < 6o — bi is equivalent 
to tpT-tiui) < '4’T-tiuo) and —Oj < bo — bi is equivalent to ipT-ti—uo) < 'ipT-ti—Ui). Since 
Uq > Ui > 0 the monotonicity of yields 

loil < 6o - bi- (30) 


An elementary calculation gives 


a'ix) 


1 

cosh(Ho + box)'^ 


i=l 


where 


fi{x) = Ui cosh(Hi + bix) cosh(Ho + box) + bi sinh(ilj + bix) cosh(ilo + box) 
— bo cosh.{Bi + bix) sinh(ilo + box). 
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The derivative of fi is 

/• (x) =bi cosh(i?o + box) (^Ui smh{Bi + bix) + (6* - bo) cosh(i?j + hx)^ 
+bo cosh.{Bi + bix) (^Ui sinh(So + box) + (bi - bo) cosh(i?o + box)^ . 
Using (30) the terms inside the brackets of each row in (31) are strictly less than 
|ai|( sgn (oj) sinh(i?j + bjx) — cosh{Bj + bjx)) < 0, {j = i, 0). 


(31) 


The last inequality is true since cosh(x) ± sinh(x) > 0. The terms outside of the brackets in 
(31) are all positive. Hence f > 0 and the fi are monotonically decreasing. Using (30) a 
simple calculation shows that hma;^_oo fi{x) = oo and lim^^-^oo fi{x) = —oo. Since > 0 for 
all i the same is true for Cie“'^/i(x). Hence g has a single maximum and is decreasing 
to the left and right of it. Furthermore g{x) > 0 and again using (30) lima;_>.oo g{x) = 
lim3;_^_oo g{x) = 0 . □ 

Proof of Lemma 4: / is continuous with compact support. Hence the extended Fourier 
transform f{z) = f {x)e~'^^ dx exists for all z G C and is analytic. For z 0,z ^ —ivk, it 
is given by 

/■K2 1 /■K2 

f{z) = / e-'^"/(x) dx = - e-'^^f{x) dx 

— il CkVk I^{vk-iz)K2 _ g(Dfe-w)Kl\ 

iz ^ Vk — \z\ ) 

k 


Since f{u — \R) = 0{u ^) for fixed R, it is absolutely integrable. By Fourier inversion 
fix) = ^ I e'^^fiz) dz = :^ r Re (e('“+«)"/(u - xR)) du, 


where the last equation follows from the fact that / is real valued and the symmetry /(z) = 
fi—z). Since f E \fiz) \ dz = JAxt\Xsi^) f 1/(2^)! dz is bounded if R G V n M, 

conditional expectation and integration can be interchanged. □ 
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